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ABSTRACT 

We report the discovery of HAT-P-14b, a fairly massive transiting extrasolar planet orbiting the 
moderately bright star GSC 3086-00152 (V = 9.98), with a period of P = 4.627669 ± 0.000005d. The 
transit is close to grazing (impact parameter 0.89lio'Qos) ano - nas a duration of 0.0912±0.0017d, with 
a reference epoch of mid transit of T c = 2,454,875. 28938±0.00047 (BJD). The orbit is slightly eccentric 
(e = 0.107±0.013), and the orientation is such that occultations are unlikely to occur. The host star is 
a slightly evolved mid-F dwarf with a mass of 1.386±0.045M Q , a radius of 1.468 ± 0.054i? Q , effective 
temperature 6600±90K, and a slightly metal-rich composition corresponding to [Fc/H] = +0.11±0.08. 
The planet has a mass of 2.232 ± 0.059 Mj and a radius of 1.150±0.052i?j, implying a mean density of 
1.82±0.24g cm -3 . Its radius is well reproduced by theoretical models for the 1.3 Gyr age of the system 
if the planet has a heavy-element fraction of about 50 (7% of its total mass). The brightness, 
near-grazing orientation, and other properties of HAT-P-14 make it a favorable transiting system to 
look for changes in the orbital elements or transit timing variations induced by a possible second 
planet, and also to place meaningful constraints on the presence of sub-Earth mass or Earth mass 
exomoons, by monitoring it for transit duration variations. 

Subject headings: planetary systems — stars: individual (HAT-P-14, GSC 3086-00152) techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

More than five dozen transiting extrasolar planets 
(TEPs) have been discovered to date, by both ground- 
based surveys and, recently, space- based surveys such 
as the CoRoT and Kepler missions (|Baglin et al.ll200& 
iBorucki et al.l l2010). The ground-based surveys are typi- 
cally able to detect TEPs orbiting only the brighter stars. 
However this disadvantage has the compensatory advan- 
tage that those TEPs detected are amenable to a wide 
range of follow-up studies, which can provide valuable 
insight into their atmospheric properties and other phys- 
ical conditions. 

Among the ground-based surveys, the Hungarian- 
made Automate d Telescope Network (HATNet; 
iBakos et al.l [200l survey has been one of the main 
contributors to the discovery of TEPs. In operation 
since 2003, it has now covered approximately 11% 
of the sky, searching for TEPs around bright stars 
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(8 < / < 12.5). HATNet operates six wide-field instru- 
ments: four at the Fred Lawrence Whipple Observatory 
(FLWO) in Arizona, and two on the roof of the hangar 
servicing the Smithsonian Astrophysical Observatory's 
Submillimeter Array, in Hawaii. Since 2006, HATNet 
has announced and published 13 TEPs. In this work we 
report our fourteenth discovery, around the relatively 
bright (V = 9.98) star we refer to as HAT-P-14, also 
known as GSC 3086-00152. 

The layout of the paper is as follows. In Section [5] we 
report the detection of the photometric signal and the 
follow-up spectroscopic and photometric observations of 
HAT-P-14. In Section[3]we examine the spectroscopic 
evidence to confirm the planetary nature of the object. 
Section @] describes the analysis of the data, beginning 
with the determination of the stellar parameters, and 
continuing with a description of our global modeling of 
the photometry and radial velocities. Results are pre- 
sented in Section [5j and discussed in Section [5] 

2. OBSERVATIONS 

2.1. Photometric detection 

The transits of HAT-P-14b were detected with the 
HAT-10 telescope located in Arizona. Observations in 
a field containing the star GSC 3086-00152 were made 
on a nightly basis between 2005 May and July, as per- 
mitted by weather conditions. We gathered 1246 expo- 
sures of 150 s with a cadence of 3.2 minutes, 206 expo- 
sures of 180 s at 3.7 minute cadence, and 1991 exposures 
of 300 s at a 5.7 minute cadence. The reason for the 
somewhat inhomogeneous data set is that the observa- 
tions were acquired during the commissioning phase of 
HAT-10. The first exposures with 150 s cadence were 
taken with an Apogee ALTA E42 back-illuminated CCD 
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Fig. 1. — Top; Unbinned /-band light curves of HAT-P-14 ob- 
tained with the HAT- 5 and HAT-10 nodes of HATNet, showing all 
3189 and 3443 instrumental measurements from fields G195 and 
G194, respectively. The photometry is phase-folded with the pe- 
riod P = 4.627669 days resulting from the global fit described in 
Section [4] The soli d lines show the transit model fit to the light 
curves (Section 14.21 1. Bottom: The same, for phases within 0.1 of 
the transit center. 

through a Cousins /-band filter. Exposure times were 
increased to 180 s after some initial testing. Finally, the 
CCD was replaced with an Apogee E10 front-illuminated 
model, and exposure times were further increased to 300 
seconds. Each image contained approximately 22,000 
stars down to a brightness limit of / ~ 14.0. For the 
brightest stars in the field we achieved a photometric 
precision of 3.5 mmag per image in the combined light 
curve. 

The calibration of the HATNet frames was carried 
out using standard photometric procedures. The cal- 
ibrated images were then subjected to star detection 
and an astrometric s olutio n, as described in more de- 
tail bv iPal &~B akos (2006). Aperture photometry was 
performed on each image at the stellar centroids de- 
rived fro m the Two Micron All Sky Survey catalog 
(2MASS; iSkrutskie etail I2006D and the individual as- 
trometric solutions. The resulting light curves were 
decorrelated (cleaned of trends) using the External Pa- 
rameter De correlation techniq ue in "constant" mode 
(EPD; see iBakos et all 120091) and the Trend Filter- 
ing Algorithm (TFA; see IKovacs et al.l 12005ft . The 
light curves were then searched for periodic box-shaped 
signals using the Box Least-Squares method (BLS; 
IKovacs et al.l 12002). We detected a significant signal 
in the light curve of CSC 3086-00152 (also known as 
2MASS 17202788+3814317; a = 17 h 20 m 27.96 s 6 = 
+38°14'31.8"; J2000; V = 9.98; IDroeee et all l2006h . 
with an apparent depth of 4.5 mmag, and a period of 
P = 4.6277 days (see Figure Q}. The drop in brightness 
had a first-to-last-contact duration q, relative the total 
period P, of q = 0.0197±0.0004, corresponding to a total 
duration of Pq = 2.189 ± 0.040 hr. 

HAT-P-14 falls in an overlapping area between two 
HATNet fields internally labeled G194 and G195. While 
the discovery was made on the basis of the G194 light 
curve, the signal was later recovered in the 3189 /-band 
observations from field G195. These latter data were 
taken during the commissioning phase of the HAT-5 tele- 
scope at FLWO, in early 2003, making this one of the 
TEPs with the longest photometric coverage. Curiously, 



the signal is present in the early HAT-5 (G195) data with 
higher significance than in field G194. The discovery of 
HAT-P-14 in 2003 was probably missed due to a com- 
bination of factors, including our different criteria for 
scheduling follow-up observations at the time, and the 
lack of the Trend Filtering Algorithm, which was only 
implemented and brought into regular use in 2005. 

2.2. Reconnaissance spectroscopy 

As is routine in the HATNet project, all candidates 
are subjected to careful scrutiny before investing valu- 
able time on large telescopes. This includes spectro- 
scopic observations at relatively modest facilities to es- 
tablish whether the transit-like feature in the light curve 
of a candidate might be due to astrophysical phenom- 
ena other than a planet transiting a star. Many of these 
false positives are associated with large radial-velocity 
variations in the star (tens of kms -1 ) that are eas- 
ily recognized. For this we have made use here of the 
Harvard- Smithsonian C enter for Astro physics (CfA) Dig- 
ital Speedometer (DS; iLathaml Il992h . an echelle spec- 
trograph mounted on the FLWO 1.5 m telescope. This 
instrument delivers high- resolution spectra (A/ A A w 
35,000) over a single order centered on the Mg lb triplet 
(~5187A), with typically low signal-to-noise (S/N) ra- 
tios that are nevertheless sufficient to derive radial veloc- 
ities (RVs) with moderate precisions of 0.5-1.0 kms -1 for 
slowly rotating stars. The same spectra can be used to 
estimate the effective temperature, surface gravity, and 
pr ojected rotati o nal ve locity of the host star, as described 
bv lTorres et al.l (|2002f) . With this facility we are able to 
reject many types of false positives, such as F dwarfs or- 
bited by M dwarfs, grazing eclipsing binaries, or triple or 
quadruple star systems. Additional tests are performed 
with other spectroscopic material described in the next 
section. 

For HAT-P-14 we obtained four observations with the 
DS between April and May of 2008. The velocity mea- 
surements showed an rms residual of 0.28kms , con- 
sistent with no detectable RV variation within the pre- 
cision of the measurements. All spectra were single- 
lined, i.e., there is no evidence for additional stars in 
the system. The atmospheric parameters we infer from 
these observations are the following: effective temper- 
ature = 6500 ± 100 if, surface gravity logg* = 
3.5 ± 0.25 (log cgs), and projected rotational velocity 
vsini = 11.7 ± 1.0 kms" 1 . The effective temperature 
corresponds to a mid-F dwarf. The mean heliocentric 
RV of HAT-P-14 is 7R V = -20.81 ± 0.28 kms" 1 . 

2.3. High resolution, high S/N spectroscopy 

Given the statistically significant transit detection by 
HATNet, and the encouraging DS results that rule out 
obvious false positives, we proceeded with the follow-up 
of this candidate by obtaining high-resolution, high-S /N 
spectra in order to characterize the RV variations with 
higher precision, and to refine the determination of the 
stellar parameters . For this we used the HIRES instru- 
ment (jVogt et al.l ll994) on the Keck I telescope located 
on Mauna Kea, Hawaii, between 2008 May 15 and 2009 
October 1. The width of the spectrometer slit was 0'.'86, 
resulting in a resolving power of A/AA « 55,000, with a 
wavelength coverage of -3800-8000 A. 
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TABLE 1 

Relative radial velocities, bisector spans, and activity 
index measurements of hat-p-14. 
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a The zero-point of these velocities is arb itrar y. An overall off- 
set 7 rc i fitted to these velocities in Section 14.21 has not been sub- 
tracted. 

b Internal errors excluding the component of astrophysical jitter 
considered in Section [5] 

c Rela tive chromospheric activity i ndex, not calibrated to the 
scale of Vaughan, Preston & Wilson (1978). 

d This observation corresponds to the iodine-free template, which 
does not yield a RV measurement but does permit the BS and S 
index to be computed. 

We obtained a total of 14 exposures through an iodine 
gas absorption cell, which was used to superimpose a 
dense forest of I2 lines on the stellar spe ctrum and estab- 
lish a n accurate wavelength fiducial (see lMarcv fc Butlerl 
1992). An additional exposure was taken without the io- 
dine cell, for use as a template in the reductions. Relative 
RVs in the solar system bary c entric frame were derived 
as described by iButler et al.l HT)9(i: , incorporating full 
modeling of the spatial and temporal variations of the 
instrumental profile. The RV measurements and their 
uncertainties are listed in Table [TJ The period-folded 
data, along with our best fit described below in Section 01 
are displayed in Figure [2J 

In the same figure we show also the relative S index, 
which is a measure of the chromospheric activity of the 
star derived from the flux in the cores of the Ca II H and 
K lines. This i ndex was computed following the pr escrip- 
tion given by iVaughan. Preston fc Wilsonl (|1978l ) , after 
matching each spectrum to a reference spectrum using 
a transformation that includes a wavelength shift and a 
flux scaling that is a polynomial as a function of wave- 
length. The transformation was determined on regions of 
the spectra that are not used in computing this indicator. 
Note that our relative S index has not been calib rated 
to the scale of IVaughan. Preston fc Wilson! ()1978[ ). We 
do not detect any significant variation of the index cor- 
related with orbital phase; such a correlation might have 
indicated that the RV variations could be due to stel- 
lar activity, casting doubt on the planetary nature of the 
candidate. There is no sign of emission in the cores of the 
Ca II H and K lines in any of our spectra, from which we 
conclude that the chromospheric activity level in HAT- 
P-14 is very low. 
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Fig. 2. — Top panel: Keck/HIRES RV measurements for 
HAT-P-14 shown as a function of orbital phase, along with our 
best-fit model (see Table [4j . Zero phase corresponds to the time 
of mid-transit. The center-of-mass velocity has been subtracted. 
Second panel: Velocity O—C residuals from the best fit. The error 
bars (which are too small to be seen in the top panel) include a 
component from astrophysical jitter (7.3ms -1 ) added in quadra- 
ture to the formal errors (see Section 0. Third panel: Bisector 
spans (BS), with the mean value subtracted. The measurement 
from the template spectrum is included. Bottom panel: Relative 
chromospheric activity index S measured from the Keck spectra. 
Note the different vertical scales of the panels. Observations shown 
twice are represented with open symbols. 

In order to permit a more accurate modeling of the 
light curve, we conducted additional photometric ob- 
servations with the KeplerCam CCD camera on the 
FLWO 1.2 m telescope. We observed five transit events 
of HAT-P-14 on the nights of 2008 October 19, 2009 
February 2, 2009 March 25, 2009 April 8, and 2009 May 
15 (Figure [3]). A Sloan i-band filter was used for all ob- 
servations. The number of images we acquired on each of 
these events is 326, 239, 470, 359, and 304, respectively. 
The exposure times were 10, 13, 20, 15, and 15 s, result- 
ing in a cadence of 22, 25, 32, 27, and 27s. The reduction 
of these images, including basic calibration, astrometry, 
aperture photo metry, and trend re moval, was performed 
as described by B akos et al.l (|2009f) . The final time series 
are shown in the top portion of Figure [3l along with our 
best-fit transit light curve model described below; the 
individual measurements are reported in Table [5] 



2.4. Photometric follow-up observations 



3. FALSE POSITIVE REJECTION 
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Fig. 3. — Unbinned instrumental transit light curves in the Sloan 
i band, acquired with KeplerCam on the FLWO 1.2 m telescope. 
The light cu rves have been processed with EPD and TFA, as de- 
scribed in S 14.21 The dates of the events are indicated. Curves 
after the first are displaced vertically for clar ity. Our best fit from 
the global modeling described in Section 14.21 is shown by the solid 
lines. Residuals from the fits are displayed at the bottom, in the 
same order as the top curves. The error bars represent the photon 
and background shot noise, plus the readout noise. 

Our initial spectroscopic analyses discussed in 
Section |2"T21 and Section |2~51 rule out the most obvious as- 
trophysical false positive scenarios. However, more sub- 
tle phenomena such as blends (contamination by an un- 
resolved eclipsing binary, whether in the background or 
associated with the target) can still mimic both the pho- 
tometric and spect roscopic signatures we see. Following 
iTorres et al.1 (j2007t ). we explore here the possibility that 
the measured Keck/HIRES radial velocities are not real, 
but are instead caused by distortions in the spectral line 
profiles due to this contamination. We determined the 
spectral line bisectors and corresponding bisector spans 
(a measure of the line asym metry) from each of our Keck 
spectra as described in §5 of lBakos et al.l ()2007l ). The bi- 
sector spans are displayed in the third panel of Figure [H 
At first glance there is little variation as a function of or- 
bital phase, in contrast to the ~400ms~ 1 peak-to-peak 



TABLE 2 

High-precision differential photometry of 
HAT-P-14 



BJD Ai a cr Ai Orig. Mag b 

(2,454,000+) (mag) (mag) (mag) 



759.57643 0.00505 0.00076 8.54130 

759.57670 0.00281 0.00076 8.53918 

759.57697 0.00561 0.00075 8.54200 

759.57725 0.00948 0.00076 8.54627 

759.57754 0.00333 0.00075 8.53897 



Note. — This table is available in a machine-readable 



form in the online journal. A portion is shown here for 

guidance regarding its form and content. 

a The out-of-transit level has been subtracted. These 

magnitudes have been subjected to the EPD and TFA 

procedures, carried out simultaneously with the transit 

fit. 

b Raw magnitude values without application of the EPD 
and TFA procedures. 

change in the RVs shown in the top panel. However, 
closer examination gives the visual impression of a slight 
correlation with the radial velocities, in the sense that 
the bisector spans for the first half of the orbit tend to 
be zero or negative (as are the RVs) , whereas those in the 
second half tend to be zero or positive (again following 
the RVs). 

15 | 1 1 1 1 1 1 1 1 1 1 
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Fig. 4. — Bisector spans ( BS) shown as a functi on of the sky con- 
tamination factor (SCF; sec Hartman et al. 2009), which is a mea- 
sure of the influence of scattered moonlight in each of our spectra 
(see text). One outlier has been removed, corresponding to the last 
Keck/HIRES observation listed in Tablc[l] which was obtained un- 
der cloudy conditions. The correlation is statistically significant, as 
indicated by a Spearman rank-order correlation false-alarm prob- 
ability of 0.04%. 

A Spearman rank-order correlation test indicates that 
a correlation as strong as that suggested between the bi- 
sector spans and the RVs would occur by chance 9.5% of 
the time (false-alarm probability, FAP), which we do not 
consider to be small enough to confirm the apparent cor- 
relation at a significant level. If we remove one outlier 
among the bisector spans near phase 0.0, correspond- 
ing to a spectrum taken through considerable clouds, 
the FAP is reduced to 6.4%, still not very compelling. 
We examined this further by considering the possibil- 
ity that the slight apparent changes in the line profiles 
might be caused by contamination not from an unre- 



HAT-P-14b 



5 



solved eclipsing binary, but from scattered moonlight un- 
der partly cloudy observing conditions. Evidence of this 
effect was seen in one of ou r previous discoveries, HAT- 
P-12 (jHartman et al J 12009( 1 . Following the same proce- 
dures used there, we computed the sky contamination 
factor (SCF) for each of our spectra of HAT-P-14, which 
quantifies the sky brightness due to moonlight and also 
accounts for the velocity difference between the moon 
and the star. We then contrasted the SCF values against 
the corresponding bisector spans. This is shown in Fig- 
ured] A rather clear trend seems to emerge, suggesting 
that the distortions in the line shapes as measured by 
the bisector spans could indeed be due to contamination 
from moonlight. A Spearman rank-order correlation test 
excluding the same outlier as before confirms this, giving 
a FAP of only 0.04%. The apparent correlation between 
the bisector spans and the velocities in Figure [2] reflects 
the fact that the six spectra taken during the second half 
of the orbit happen to have SCF values that are larger, on 
average, than those in the first half. Based on the above, 
we conclude that the small changes in the line shapes as 
measured by the highly sensitive bisector spans are well 
explained by contamination from scattered moonlight, 
and that there is no real correlation with the radial ve- 
locities. The radial velocity variations in HAT-P-14 are 
therefore real, and must be due to a planetary companion 
in orbit around the star. 

4. ANALYSIS 
4.1. Properties of the parent star 

Fundamental parameters of the host star HAT-P-14 
such as the mass (M*) and radius (R*), which are needed 
to infer the planetary properties, depend strongly on 
other stellar quantities that can be derived spectroscop- 
ically. For this we have relied on our template spec- 
trum obtained with the Keck/HIRES instrument, and 
the an alysis package known as S pectroscopy Made Easy 
(SME; iValenti & Piskunovlll996l) a long w ith the atomic 
line database of IValenti fc Fischer! ([2005) . This analysis 
yielded the following initial values: effective tempera- 
ture T e ff* = 6310 ± 80 K, surface gravity log = 3.80 ± 
0.10 (cgs), metallicity [Fe/H] = -0.02±0.05dex, and pro- 
jected rotational velocity vsini = 8.8 ± 0.5 krns" 1 . The 
temperature and metallicity uncertainties have been con- 
servatively increased by a factor of two over their nomi- 
nal values, based on previous experience, to account for 
possible systematic errors. 

In principle the effective temperature and metallicity, 
along with the surface gravity taken as a luminosity in- 
dicator, could be used as constraints to infer the stellar 
mass and radius by comparison with stellar evolution 
models. However, the effect of log g+ on the spectral line 
shapes is rather subtle, and as a result it is typically dif- 
ficult to determine accurately, so that it is a rather poor 
luminosity indicator for this particular application. A 
trigonometric parallax is unfortunately unavailable for 
HAT-P-14 since the star was n ot included among th e 
targets of the Hipparcos mission (jPerrvman et al1ll997(l . 
despite being bright enough. In any case, for planetary 
transits a stronger constraint on the luminosity is of- 
ten provided by a/i?*, the normalized semimajor axis, 
which is closely related to p i; the mean stellar den- 
sity. The quantity a/i?* can be derived directly from 



TABLE 3 
Stellar parameters for HAT-P-14 



Parameter Value Source 



Spectroscopic properties 

(K) 6600 ± 90 SME a 

[Fe/H] +0.11 ± 0.08 SME 

vsini (kms- 1 ) . . . 8.4 ±0.5 SME 

Umac (kms" 1 ) .... 5.30 SME 

v mic (kms" 1 ) .... 0.85 SME 

Trv (kms" 1 ) -20.81±0.28 DS 

Photometric properties 

B T (mag) 10.494 ± 0.031 Tycho-2 

V T (mag) 10.038 ± 0.029 Tycho-2 

V (mag) 9.980 ± 0.058 TASS 

I c (mag) 9.453 ± 0.102 TASS 

J (mag) 9.094 ± 0.021 2MASS 

H (mag) 8.927 ± 0.020 2MASS 

K„ (mag) 8.851 ± 0.019 2MASS 

E(B — V) (mag) . . 0.028 ± 0.015 Dust maps b 

Derived properties 

Af* (M @ ) 1.386 ±0.045 YY+a/fl*+SME c 

R* (Rq) 1.468 ±0.054 YY±a/.R*±SME 

log g+ (cgs) 4.25 ±0.03 YY±a/R*±SME 

L* (L Q ) 3.66 ±0.37 YY+a/ii*+SME 

M v (mag) 3.31 ±0.12 YY±a/i?*±SME 

Age (Gyr) 1.3 ±0.4 YY±a/i?*±SME 

Distance (pc) 205 ± 11 YY±a/i?*±SME 



a SME — "Spectroscopy Made Easy" package for the analysis of high- 
resolution spectra (Valcnti & Piskunov 1996). These parameters rely 
primarily on SME, but have a small dependence also on the iterative 
analysis incorporating the isochronc search and global modeling of the 
data, as described in the text. 

Average of results from ISchlegel et al.l dl998f) . IBurstein &: Heilesj 
J 1982fl . and lHakkila et alJ fl997D . with small corrections for distance. 
c YY+o/fl,+SME = Based on the YY isochrones flYi et al.ll200ll) , 
a/R* as a luminosity indicator, and the SME results. 

the tra nsit l ight curves (see lSozzetti et al.1l2007l and also 
Section |4"T2"|) . This, in turn, allows us to improve on the 
determination of the spectroscopic parameters by sup- 
plying an indirect constraint on the weakly determined 
spectroscopic value of log g+ that removes degeneracies. 
We take this approach here, as described below. 

Our initial values of T e g+, log 3*, and [Fe/H] were used 
to determine auxiliary quantities needed in the global 
modeling of the follow-up photometry and radial veloc- 
ities (specifically, the limb-darkening coefficients). This 
modeling, the details of which are described in the next 
section, uses a Monte Carlo approach to deliver the nu- 
merical probability distribution of a/R+ and other fitted 
variable s. For further details we refer the reader to iPall 
(2009b). When combining a/i?* (used as a proxy for lu- 
minosity) with assumed Gaussian distributions for T c g* 
and [Fe/H] based on the SME determinations, a compar- 
ison with stellar evolution models allows the probability 
distributions of other stellar properties to be inferred, 
including logg*. Here we use the stellar evol ution cal- 
culati ons from the Yonsei-Yale (YY) series by lYi et al.l 
(20031). The comparison against the model isochrones 
was carried out for each of 10,000 Monte Carlo trial sets 
(see Section FO]) . Parameter combinations correspond- 
ing to unphysical locations in the H-R diagram (1% of 
the trials) were ignored, and replaced with another ran- 
domly drawn parameter set. The result for the surface 
gravity, log g+ = 4.25 ± 0.03, is significantly different 
from our initial SME analysis, which is not surprising 
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in view of the strong correlations among T ff*, [Fe/H], 
and log g+ that are often present in spectroscopic de- 
terminations. Therefore, we carried out a second iter- 
ation in which we adopted this value of log g+ and held it 
fixed in a new SME analysis (coupled with a new global 
modeling of the RV and light curves), adjusting only 
T off *, [Fe/H], and vsxai. This gave T efr * = 6600 ± 90K, 
[Fe/H] = +0.11 ± 0.08, and vsini = 8.4 ± O-Skms" 1 , 
in which the conservative uncertainties for the first two 
have been increased by a factor of two over their for- 
mal values, as before. A further iteration did not change 
log significantly, so wc adopted the values stated above 
as the final atmospheric properties of the star. They are 
collected in Table [3J together with the adopted values for 
the macroturbulent and microturbulent velocities. 

We note that the effective temperature also changed 
rather significantly from the first SME iteration, by 
about 300 K. An external check on the accuracy of T e ff* 
may be obtained from color indices derived from ex- 
isting brightness measurements for HAT-P-14. Photo- 
metric measurements on a standard system are avail- 
able from the T y cho-2, 2MASS and T ASS catalogs 
(lH0g et all 120001: iSkrutskie etaLl 120061: iDroege et all 
2006), and are summarized in the second block of 
Table [3J We constructed seven different color in- 
dices, and applied color /t emperature calibrations by 
IRamrrez fc Me lcndcz (2005) (avail able for six of the color 
indi ces). iCasagrande et al.l (I2006D (four col or indices), 
and iGonzalez H ernandez & Bonifacio (2009) (five color 
indices). These relations include terms that depend on 
metallicity. Because such photometric temperatures are 
ver y sensitive to r e ddeni ng, we examined the dust maps 
by ISchlegel et all (|1998D in the direction to HAT-P-f 4 
and obtained E{B — V) — 0.030, after a small correction 
for distance using the value of 205±11 pc reported below. 
Two other reddening esti mates of 0.022 and 0.033 were 
obtained simi l arly from IBurstein fc Heilell ([1982D and 
IHakkila et al.l (|1997t ). respectively. We adopt the straight 
average of the three measures, E(B—V) = 0.028 ± 0.015, 
with a conservative uncertainty. 

We corrected each of the color indices for reddening, 
and computed weighted temperature averages for each 
set of calibrations, with weights determined from the in- 
dividual temperature uncertainties. The latter include 
all photometric errors, the uncertainty in both the red- 
dening and the metallicity, and the scatter of the cal- 
ibrations themselves. The seven color indices are not 
completely independent, but nevertheless give a useful 
sense of the scatter. The agreement between the 15 sep- 
arate estimates from all three calibrations is quite sat- 
isfactory, the dispersion being 112K. The mean of the 
three calibration averages is T e g = 6598 ± 100 K, which 
is in excellent accord with the final spectroscopic value 
from SME and supports its accuracy. 

With the adopted spectroscopic parameters the model 
isochrones yield the stellar mass and radius M* = 1.386± 
0.045 M Q and i?* = 1.468 ± O.O54i? , along with other 
properties listed at the bottom of Table [3J HAT-P-14 
is a slightly evolved F star with an estimated age of 
1.3 ± 0.4 Gyr, according to these models. The inferred 
location of the star in a diagram of a/i?+ versus T^g*, 
analogous to the classical H-R diagram, is shown in Fig- 
ure [5] The stellar properties and their ler and 2<r con- 



fidence e l lipsoi ds are displayed against the backdrop of 
lYi et al.l (|2001l) isochrones for the measured metallicity 
of [Fe/H] = +0.11, and a range of ages. For comparison, 
the location implied by the initial SME results is also 
shown (triangle), and corresponds to a somewhat more 
evolved state. 




7000 6500 6000 5500 5000 



Effective temperature (K) 

Fig. 5. — Model isochrones from Yi et al. (2001) for the measured 
metallicity of HAT-P-14, [Fe/H]= +0.11, and ages of 0.2, 0.4, 0.6, 
0.8, 1.0, 1.5, 2.0, 2.5, and 3.0 Gyr (left to right). The adopted 
values of T c g+ and a/R* are shown together with their la and 2cr 
confidence ellipsoids. The initial values of T^ff* and a/R+ from the 
first SME and light curve analyses are represented with a triangle. 



The distance to the star may be obtained by compar- 
ing the absolute magnitudes inferred from the stellar evo- 
lution models against the apparent brightness measure- 
ments mentioned earlier. For this we use the V, Ic, 
J, H, and K s magnitudes, corrected for extinction using 
A(V) = 3.1 x E(B—V) appropriately scaled to each pass- 
band following |Cardi!iretaL| $389). The near-infrared 
magnitudes from 2MASS were transformed t o the ESQ 
system of the isochrones using the relations bv lCarpenterl 
([20011 ) . The five distance estimates agree with each other 
to within a few pc. The average is 205 ± 11 pc, where 
the uncertainty includes all sources of error (photomet- 
ric and spectroscopic) , but excludes possible systematics 
from the evolution models themselves, which are difficult 
to quantify. 

4.2. Global modeling of the data 

This section describes the procedure we followed to 
model the HATNet photometry, the follow-up photome- 
try, and the radial velocities simultaneously. Our model 
for t he follow-up light cur ves used analytic formulae 
from iMandel fc Agoll ([2002) for the eclipse of a star by 
a planet, with limb darkening being prescribed by a 
quadratic law. The limb darkening coefficients for the 
Sloan i band were interpolated from the tables bv lClaretl 
(2004) for the spectroscopic parameters of the star as 
determined from the SME analysis (Section [OJ . The 
transit shape was parametrized by the normalized plan- 
etary radius p = R p /R*, the square of the impact pa- 
rameter b 2 , and the reciprocal of the half duration of 
the transit (/R±. We chose these parameters because of 
their simple geometric meanings and the fact tha t they 
show negligible correlations (see lBakos et al. 2009). The 
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relation between C/R* and the quantity a/R+, used in 
Section [4~T1 is given by 

a/R* = P/2n{(/R+)^l - b 2 y/l - e 2 /(l + esinw) (1) 

(see, e.g., iTinglev fc Sac kctt 20Q5|). Our model for the 
HATNet data was the simplified "P1P3" version of the 
iMandel fc Agoll ()2002|) analytic functions (an expansion 
in terms Legendre pol y nomia ls), for the reasons de- 
scribed by IBakos et al.l (120091) . Following the formal- 
ism presented bv lPall (|2009D . the RVs were fitted with 
an eccentric Keplerian model parametrized by the semi- 
amplitude K and Lagrangian elements k = e cos u and 
h = esinw, in which ui is the longitude of periastron. 

We assumed that there is a strict periodicity in the 
individual transit times. We assigned the transit num- 
ber Ntr — to the first complete follow-up light curve 
gathered on 2009 March 25. The adjustable parameters 
in the fit that determine the ephemeris were chosen to 
be the time of the first transit center in HATNet field 
G195, T Ci _473, and that of the last transit center observed 
with the FLWO 1.2m telescope, T Cj+ n. We used these 
as opposed to the period and reference epoch i n order to 
minim ize correlations between parameters (see iPal et ahl 
2008). Times of mid-transit for intermediate events were 
interpolated using these two epochs and the correspond- 
ing transit number of each event, N tr . The eight main pa- 
rameters describing the physical model were thus T c ._473, 
T c + n, R p /R+, b 2 , C/ii*, K, k = ecosuj, and h = esinw. 
Five additional ones were included that have to do with 
the instrumental configuration. These are the HATNet 
blend factors -Bi ns t (one for each HATNet field, G194 and 
G195), which account for possible dilution of the transit 
in the HATNet light curves from background stars due to 
the broad PSF (20" FWHM), the HATNet out-of-transit 
magnitudes iW^HATNet ( one f° r each field), and the rela- 
tive zero-point 7 rc i of the Keck RVs. 

We extended our physical model with an instrumen- 
tal model that describes brightness variations caused 
by systematic errors in the measurements. This was 
done in a similar fashion to the analysis presented by 
IBakos et all (|2009h . The HATNet photometry had al- 
ready been EPD- and TFA-corrected before the global 
modeling, so we only considered corrections for system- 
atics in the follow-up light curves. We chose the "ELTG" 
method, i.e., EPD was performed in "local" mode with 
EPD coefficients defined for each night, and TEA was 
performed in "global" mode using the same set of stars 
and TFA coefficients for all nights. The five EPD pa- 
rameters were the hour angle (representing a monotonic 
trend that changes linearly over time), the square of the 
hour angle (reflecting elevation), and the stellar profile 
parameters (equivalent to the FWHM, elongation, and 
position angle of the image). The functional forms of the 
above effects contained six coefficients, including the aux- 
iliary out-of-transit magnitude of the individual events. 
The EPD parameters were independent for all five nights, 
implying 30 additional coefficients in the global fit. For 
the global TFA analysis we chose 20 template stars that 
had good quality measurements for all nights and on 
all frames, implying an additional 20 parameters in the 
fit. Thus, the total number of fitted parameters was 13 
(physical model with 5 configuration-related parameters) 
+ 30 (local EPD) + 20 (global TFA) = 63, i.e., much 



smaller than the total number of data point s (1711). 

The joint fit was performed as described in lBakos et al.l 
(2009). We minimized x 2 m the space of parame- 
ters by using a hybrid algorithm, co mbining the down- 
hill simplex method (AMOEBA; see iPress et al.l Il992t ) 
with a classical least squares algorithm. Uncertainties 
for the parameters were derived applyi ng the Ma rkov 
Chain Monte-Carlo method CM CMC. see lFordll2006fl us- 
ing "Hyperplane-CLLS" chains (jBakos et al.ll2009ft . This 
provided the full a posteriori probability distributions 
of all adjusted variables. The length of the chains was 
10,000 points. The a priori distributions of the param- 
eters for these chains were chosen to be Gaussian, with 
eigenvalues and eigenvectors derived from the Fisher co- 
variance matrix for the best-fit solution. The Fisher 
covariance matrix was calc ulated ana lytically using the 

partial derivatives given bvipm (pool . 

5. RESULTS 

With the procedure above we obtained the a posteriori 
distributions for all fitted variables, and other quantities 
of interest such as a/fi t . As described in Section |4~T1 
a/R± was used together with stellar evolution models to 
infer a theoretical value for logg* that is significantly 
more accurate than the spectroscopic value. The im- 
proved estimate was in turn applied to a second iter- 
ation of the SME analysis, as explained previously, in 
order to obtain better estimates of T c s* and [Fe/H]. The 
global modeling was then repeated with updated limb- 
darkening coefficients based on those new spectroscopic 
determinations. The resulting geometric parameters per- 
taining to the light curves and velocity curves are listed 
in Table [4] The rms residual of the RV measurements 
from the orbital fit is 8.1ms _1 . 

Included in this table is the RV "jitter". This is a 
component of assumed astrophysical noise intrinsic to the 
star that we added in quadrature to the internal errors 
for the RVs in order to achieve x 2 /dof = 1 from the RV 
data for the global fit. Auxiliary parameters not listed in 
the table are: T c _ 473 = 2,452,728.0512 ± 0.0024 (BJD), 
T c +ii = 2,454,967.8427 ± 0.0005 (BJD), the blending 
factors Binstr = 0.73 ± 0.09 for HATNet field G194 and 
B instI = 0.96 ± 0.04 for field G195, and 7rc i = 26.6 ± 
2.4ms -1 . The latter quantity represents an arbitrary 
offset for the Keck RVs, and does not correspond to the 
true center of mass velocity of the system, which was 
listed earlier in Table [3] (trv)- 

The planetary parameters and their uncertainties can 
be derived by combining the a posteriori distributions for 
the stellar, light curve, and RV parameters. In this way 
we find a mass for the planet of M p = 2.232 ± 0.059 Mj 
and a radius of R p — 1.150 ± 0.052 leading to a mean 
density p p — 1.82 ± 0.24gcm -3 . These and other plane- 
tary parameters are listed at the bottom of Table |4j 

We note that the eccentricity of the orbit is signifi- 
cantly different from zero: e = 0.107 ± 0.013. With a 
longitude of periastron u = 95° ± 4°, the orientation of 
the system is such that the orbit is viewed very nearly 
along the major axis, with periastron on the near side. 
Furthermore, the eccentricity is large enough that the 
separation between the star and the planet when the 
planet is behind is some 23% larger than at transit. 
This, combined with the relatively low inclination angle, 
implies that there should be no occultations of HAT-P- 
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TABLE 4 
Orbital and planetary parameters 

Parameter Value 

Light curve parameters 

P (days) 4.627669 ± 0.000005 

T c (BJD) a 2,454,875.28938 ±0.00047 

Ti 4 (days) a 0.0912 ± 0.0017 

T 12 = T 34 (days) a 0.0287 ± 0.0026 

a/R+ 8.87 ±0.29 

C/-R* 29.53 ±0.42 

Rp/R* 0.0805 ± 0.0015 

0.794±°;°i* 

6 = acosi/-R* 0.89lto°og 

i (deg) 83.5 ± 0.3 

Limb-darkening coefficients (Sloan i band) b 

ci (linear term) 0.1089 

C2 (quadratic term) 0.2439 

RV parameters 

K (ms" 1 ) 219.0 ±3.3 

fc = ecosw c -0.009 ±0.009 

h = esinw c 0.106 ± 0.013 

e 0.107 ±0.013 

w (deg) 94 ± 4 

RV jitter (ms -1 ) 7.3 

o-rv (ms -1 ) 8.1 

Planetary parameters 

M p (Mj) 2.232 ±0.059 

R p (i?j) 1.150 ±0.052 

C(M p ,Rp) d 0.45 

p p (gem" 3 ) 1.82 ±0.24 

logg p (cgs) 3.621 ±0.037 

a (AU) 0.0606 ± 0.0007 

T eq (K) 1570 ± 34 

6° 0.168 ±0.008 

(F) (10 9 ergs" 1 cm~ 2 ) f 1.37 ±0.12 

a T c : Reference epoch of mid transit that minimizes the cor- 
relation with the orbital period. It corresponds to Ntr — — 9- 
T14: total transit duration, i.e., the interval of time between 
the first and last contacts; T12 — T34: ingress/egress time, 
i.e., the time interval between the first and second, or third 
and fourth contacts. 

k Values for a quadratic law, adopted from the tabulations by 
Claret ( 2004) according to the spectroscopic (SME) parame- 
ters listed in Table [3] 

c Lagrangian orbital parameters derived from the global mod- 
eling, and primarily determined by the RV data. 

Correlation coefficient between the planetary mass M p and 
radius R p . 

e Safronov number given by — ^(V esc/Vorb) 2 — 
(a/R p )(M p /Mi,) fsee lHansen &: Barrhanll2007l) . 
* Incoming flux per unit surface area, averaged over the 
slightly eccentric orbit. 



14b. Formally, the impact parameter at the occultation 
is ^occult = 1.102 ± 0.030, which is greater than unity at 
the 3<t confidence level. 

Individual times of mid transit were measured for each 
of the five events observed with KeplerCam. This was 
done by holding fixed all model parameters that define 
the shape of the transit to their final values in Table [4] 
and adjusting only the time of the center of the event 
in each series. The results are listed in Table [SJ along 
with their uncertainties and O — C residuals from the 
ephemeris for the planet. As expected, the events with 
only partial coverage of the transit (first two in Table [SJ 
see also Figure [3]) have larger uncertainties and residuals. 

As a final note, the barycentric Julian dates in this 



TABLE 5 

Transit timing measurements for HAT- P- 14 



Ntr 


(BJD-2,454,000) 


(days) 


O-C 
(sec) 


-34 


759.59280 


0.00290 


-458 


-11 


866.03268 


0.00273 


-127 


+0 


916.93940 


0.00081 


±90 


±3 


930.82160 


0.00114 


+24 


+11 


967.84186 


0.00088 


-59 



table and throughout the paper are calculated from Co- 
ordinated Universal Time (UTC), in which a second has 
the same duration as a second of International Atomic 
Time (TAI), but the scale is subject to occasional 1- 
second adjustments (leap seconds) to stay in pace with 
the slowly and irregularly decreasing rotation rate of the 
Earth. UTC (the basis of civil time) is therefore not con- 
tinuous, and is unsuitable as a physical time coordinate. 
This can be important for some applications that com- 
bine timing measureme nts over long periods of time (e.g., 
IGuinan fe" Ribas 2001), a problem that has long been 
recognized in the binary sta r community (for a pointed 
discussion, see Bastian 2000) as well as the Solar System 
and pulsar communities. Timing measurements in the 
transiting planet field now span a decade, and although 
the differential effect over this period is only a few sec- 
onds (significantly smaller than the typical measurement 
precision), it will become more important in the future. 
The issue has already been mentioned recently in the con- 
text of s tudying transit tim ing variations for transiting 
planets (jAdams et al.l 12010) . While the use of Terres- 
trial Time (TT = TAI + 32.184 seconds) as the basis for 
Julian date calculations eliminates the problem, for prac- 
tical reasons virtually all the JDs (and BJDs) reported 
in previous transiting planet studies, including those in 
the HATNet series, are implicitly on the UTC scale. We 
continue this practice here, but caution the reader that 
the TT— UTC corrections need to be taken into account 
for the most demanding applications. 

6. DISCUSSION 

At V = 9.98, HAT-P-14 is among the brightest TEP 
host stars presently known (only 11 of the other 66 
are brighter), which should facilitate a wide range of 
follow-up studies. Our planetary mass determination of 
M p = 2.23 Mj places HAT-P-14b in the category of 
the massive transiting planets discovered to date, which 
seem to be less common than those with smaller masses. 
In fact, a look at the mass distribution for the known 
transiting systems shows that it falls off precisely around 
2 Mj . Since there is no obvious selection effect that would 
make these massive objects more difficult to detect or 
confirm (such as smaller radii, longer orbital periods, 
or fainter or hotter parent stars), it may be that they 
are indeed less co mmon as a population. It has also 
been noticed (e.g., iSouthworth et al.l 120091 : Uoshi et all 
I2009f) that the most eccentric cases among the transit- 
ing planets are all massive: all systems with e > 0.2 
are more massive than 3Afj, although the numbers are 
still small. A similar trend is seen among the non- 
transiting planets, bearing in mind the sini ambiguity in 
their masses. Among the transiting systems whose spin- 
orbit alignment has been measured through observation 
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of the Rossiter-McLaughlin effect, it would also appear 
that the more massive ones are disproportionately mis- 
aligned compared to the ones with smaller-mass planets, 
considering their relative populations, although again the 
numbers are as yet too small to be conclusive. Obviously 
it would be useful to determine the spin-orbit alignment 
of the HAT-P-14 system, which should be readily mea- 
surable given that vsini — 8.4 km s -1 . 

The mass of HAT- P - 14b is quite simil ar to the recently 
discovered Kepler-5b (jKoch et al.ll2010l ) , but its radius is 
about 20 % smaller. A c ompa rison with the evolution 
models of iFortnev et al.l (|2007f ). at an equivalent solar 
semimajor axis of a equ iv = 0.0317 AU, indicates that the 
radius of HAT-P-14b is well reproduced for the 1.3 Gyr 
age of the system if it has about 50 M® worth of heavy 
elements in its interior (about 7% of its total mass). This 
amount of metals is consistent with the correlation be- 
tween co re mass and metalli city of the parent star pro - 
posed bv lGuillot etHI (|2006t) and lBurrows et al.l (|2007l ). 
which seems to support the core-accretion mode of planet 
formation. The incident flux we compute for the planet 
averaged over its eccentric orbit, (F) = 1.37 ±0.12 x 10 9 
ergs -1 cm" 2 , places it in the proposed pM category of 
IFortnev et al.l (|2008fl . These objects are expected to 
present temperature inversions in their atmospheres, and 
large day/night temperature contrasts. The equilibrium 
temperature of HAT-P-14b, assuming zero albedo and 
full redistribution of the incident radiation, is 1570±34K. 

The non-zero eccentricity of the orbit for the relatively 
short period of 4.6 days raises the question of whether 
tidal forces have had enough time to circularize the or- 
bit at the 1.3 ± 0.4 Gyr age of the system, or alterna- 
tively, whether there might be a second planet perturbing 
the orbit of HAT-P-14b and pumping up its eccentricity. 
Following lAdams fc Laughlinl (|2006), we estimate the 
timescale for tidal circularization to be r C j rc w 0.33 Gyr, 
although this is a strong function of the poorly known 
tidal quality factor Q p , for which we have adopted here 
the commonly used value of 10 5 . For Q p ~ 10 6 the 
timescale is 10 times longer, and a recent study of the 
0.78-day transiting planet WASP- 19b has sug gested Q p 
could be much larger still (|Hebb et al.l 120101 ). Conse- 
quently we cannot rule out that the eccentricity we mea- 
sure for HAT-P-14b is primordial, at least in part, nor 
can we exclude the presence of a second planet in the 
system. 

The grazing orientation (b — 0.891+o;QQg) signifi- 
cantly enhances the detection sensitivity to additional 
planets in the system. These could manifest them- 
selves through transit duration variations (TDVs), or 
equivalently, changes in the orbital inclination angle 
as has been claimed for TrES-2 b by iMislis fc Schmittl 
(I2009aft and IMislis et all (|2009bl) . and recently disputed 
bv lScuderi et al.l (|2010l ). Perturbing planets may also in- 
duce variations in the times of mid-transit (transit timing 
variat ions, or TTVs) through grav i tational interaction 
(e.g., iDobrovolskis fc Boruckil 119961 : iHolman fc Murray! 
I2005D . In general, there exists no exact closed- form ex- 
pressions for evaluating the amplitudes of such signals 
since no general solution to the 3-body problem exists. 
One typical approach i s to employ an ap proximation to 
simplify the problem. lAgol et all ((2005) provided use- 
ful analytic expressions for the case of a coplanar sys- 



tem where the host star has a much larger mass than 
the two companion planets, for several different scenar- 
ios. The largest TTVs occur for planets in mean mo- 
tion resonance, and we may use expressions [33] and [34] 
from lAgol et all ((2005) to estimate the amplitude and 
libration period of the TTV signal in a 2:1 resonance. 
For the HAT-P-14 system we find that such a planet 
would induce a peak-to-peak deviation in the timings of 
13.6s and 126.3s for a mass of 0.1 M® and 1M®, re- 
spectively, and the libration period would be ~175 days. 
Thus, an Earth-mass planet in this configuration would 
be eminently detectable with further observations. Out- 
of-reso nance perturbers c an be considered using equation 
[32] of lAgol et ail (|2005[ ). but we find they lead to sub- 
second amplitudes for various configurations, far beyond 
the current measurement capabilities. 

The properties of HAT-P-14b also make it attractive 
for a future search for a companion satellite ('exomoon'). 
The Hill radius (Rh) extends to 8.1 planetary radii 
9.1 R j), and using the expressions of iDomingos et all 
2006), which account for the planet's orbital eccentric- 
ity, the maximum distance at which an exomoon could 
be stable for this planet would be 0.436 Rh and 0.824 Rh 
for a prograde and retr ograde exomoon, respectively. 
I Barnes fc O' Brier] (|2002f ) presented analytic approxima- 
tions for the maximum stable exomoon mass for close-in 
extrasolar giant planets based upon tidal dissipation ar- 
guments, which were shown to provide excellent agree- 
ment with numerical integrations. Assuming again a 
tidal dissipation value of Q p ~ 10 5 and a Love number 
of A)2p — 0.51, the maximum stable moon mass over the 
1.3 Gyr lifetime of the star is 0.002 M® and 0.12 Af® for a 
moon at the maximum planet-moon separations for pro- 
grade and retrograde orbits, respectively. For Q p ~ 10 6 
these limits would be 10 times larger, or 1.2 M® in the 
retrograde case. Such exomoons would produce small 
TTVs, but the highly gr azing orientation provides a bet- 
ter means of detection. iKippingj (|2009al lbh showed that 
an exomoon induces two types of transit duration varia- 
tions on the host planet: a sky-projected tangential ve- 
locity variation (the V-component), and a spatially or- 
thogonal transit impact parameter variation (the TIP 
component). While the TTV and TDV-V rms ampli- 
tudes of a 1.2 M® exomoon would be 4.5 s and 1.0 s, re- 
spectively, the TDV-TIP component could be as high 
as 16 s or 45 s peak-to-peak for the optimum configu- 
ration of a favorably inclined exomoon. For compari- 
son, another near-grazi ng transiting planet is TrES-2b 
(|Q 'Donovan et al.l 12006). and repeating the same calcu- 
lation (with Q p ~ 10 6 , as above) we find that the max- 
imal TDV-TIP of the maximum-mass stable retrograde 
exomoon would be only 0.07 s. The much smaller value 
in this case is mostly due to the low maximum stable 
moon mass of 0.005 M®, as a result of the lower plane- 
tary mass and longer system age. We suggest, therefore, 
that future monitoring of the transit duration of HAT- 
P-14b could allow for interesting constraints not only on 
the presence of additional planets, but also of sub-Earth 
mass or Earth mass exomoons. 
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